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We employ our new approach to non-relativistic supersymmetric quantum mechanics (SUSY- 
QM), (J. Phys. Chem. A 114, 8202(2010)) for any number of dimensions and distinguishable 
particles, to treat the hydrogen atom in full three-dimensional detail. In contrast to the standard one- 
dimensional radial equation SUSY-QM treatment of the hydrogen atom, where the superpotential is 
a scalar, in a full three-dimensional treatment, it is a vector which applies regardless of the electron 
angular momentum. The original scalar Schrodinger Hamiltonian operator is factored into vector 
"charge" operators: Q and Q'. Using these operators, the first sector Hamiltonian is written as 
Hi = ■ Q + Eq . The second sector Hamiltonian is a tensor given 

and is isospectral with rti. The second sector ground state, tpQ , can be used to obtain the excited 
state wave functions of the first sector by application of the adjoint charge operator. Alternatively, 
Q applied to analytical, sector one excited states yield analytical results for the sector two vector 
eigenstates. Several of these are plotted for illustration. We then adapt the aufbau principle to 
show this approach can be applied to treat the helium atom. 

PACS numbers: 31.15.A-, 11.30.Pb 



I. INTRODUCTION 

In a previous publication, we have provided a gener- 
alization of supersymmetric quantum mechanics (SUSY- 
QM) to treat any number of dimensions or particles with 
a focus on its usefulness as a computational tool for calcu- 
lating accurate excited state energies and wave functions 
[T] . We note that Stedman has presented a similar treat- 
ment for multi-dimensional systems. However, he gener- 
ates extra sector states not present in our approach|2j. 
His result is more complicated and the meaning of the 
sector Hamiltonians is not discussed. Stedman presents 
his equations for the hydrogen atom but does not solve 
any of the higher sector equations. Because of the signifi- 
cant analytical and computational ramifications, we here 
apply our multi-dimensional generalization of SUSY-QM 
to the hydrogen atom in full three-dimensional detail. 
This is of interest because, until now, the standard ap- 
plication of SUSY-QM to the hydrogen atom required 
that we first separate out the angular degrees of freedom 
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- effectively reducing the problem to a one-dimensional 
treatment (3HS] . With our vector superpotential ap- 
proach, one can deal with the full three-dimensional na- 
ture of the hydrogen atom. 

Our approach provides, for the first time, a SUSY-QM 
framework that can be employed to treat non-hydrogenic 
atoms. For example, the standard SUSY-QM treatment 
of the hydrogen atom cannot be readily extended to the 
helium atom because it is impossible reduce it to a one- 
dimensional system. In addition, the form of the three- 
dimensional vector superpotential for the hydrogen atom 
is of interest in its own right. It is quite different from 
the radial superpotential obtained in earlier SUSY-QM 
studies of the hydrogen atom. The present study thus 
lays the groundwork for a systematic SUSY-QM study of 
excited state energies and wave functions of atoms. As 
we have discussed in our earlier studies, there appears to 
be a significant increase in the accuracy of both excited 
state energies and wave functions when one computes 
the sector two ground state energy and wave functions, 
followed by application of the adjoint "charge operator" 
to generate the sector one excited state wave functions 

BEIE]. 

This paper is organized as follows. In section II, to 
make this paper self contained, we give a detailed descrip- 
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tion of the SUSY-QM formulation valid for any number 
of dimensions and particles. In section III, we describe 
the application of the theory to the hydrogen atom in 
three spatial dimensions. We give exact analytical ex- 
pressions for the charge operators and obtain the ground 
state wave functions and energies of the sector two hy- 
drogen atom tensor Hamiltonian. In fact, we note that 
since the exact solutions of the sector one hydrogen atom 
Hamiltonian are analytically known, we can also gener- 
ate any of the excited states of the sector two problem. 
In Section IV, we outline how the approach can be ap- 
plied to the helium atom. In Section V, we examine the 
ramifications of the aufbau principle in the second sec- 
tor. Finally, in Section VI, we present our conclusions 
and plans for future research. 



II. INTRODUCTION TO THE PARTNER 
HAMILTONIAN FORMULATION OF SUSY-QM 
IN N-DIMENSIONS 

The standard approach to supersymmetric quantum 
mechanics (SUSY-QM) provides an elegant scheme for 
solving one-dimensional problems, but until recently, it 
has not been generalized to multiple dimensions (5MlO|. 
We were originally attracted to SUSY-QM as a novel 
computational approach, since it converts a standard 
second order differential equation to one of first order. 
This reduction of order is achieved by factoring the 
Schrodinger Hamiltonian operator in terms of so-called 
"charge" operators, Q and QK The simplest, and best 
known example is the one-dimensional harmonic oscil- 
lator where the Q and are the well-known lowering 
and raising operators [TT]. We have solved a number of 
one-dimensional problems using this approach. Because 
of the success of our one-dimensional studies, we were 
motivated to generalize SUSY-QM to any number of di- 
mensions or particles. 

Most previous attempts to generalize SUSY-QM to 
treat more than one spatial dimension and more than one 
particle generally have involved introducing additional 
"spin-like" degrees of freedom |12ff2"T] . In our method [lj, 
we make use of a vectorial approach that simultaneously 
treats more than one dimension and any number of dis- 
tinguishable particles (see also |2])- We consider, there- 
fore, a system of n-particles in three-dimensional space. 
We denote the coordinates of particle i by (xi, Ui,Zi). We 
then define an orthogonal hyperspace of dimension 3n. 
We take the Hamiltonian for this system to be given by 

-Hi = -V 2 + Ui (III) 

where 

* = (IL2) 

3 J 

and tj ■ <ffc = 8jk- The subscript "1" indicates this is 



the "sector one" Hamiltonian. For simplicity we take the 
masses of the particles to be equal and use units such 
that h 2 /2m — 1. For the development here we assume a 
Cartesian coordinate space, but have provided an exten- 
sion to more general curvilinear coordinates in a previous 
publication [1]. 

As per usual in quantum mechanics, the ground-state 
wave function is a solution of the Schrodinger equation, 

Hx^=E^\ (II.3) 

We also emphasize that the lowest energy state, \ is 
nodeless. 

We now define a vector superpotential, W, as 

W = -fln^ 1) , (IL4) 

which is to say 

3n 3n ~ 

^ = E ?i w i = - E 3 ln ^ ] ■ ( IL5 ) 

3=1 3=1 3 

It is straightforward to see that one can write Hi in terms 
of W as 

= (-d i + w i )(d i + w i ), 

where, according to the Einstein convention, we sum over 
repeated indices. Since (V + W)^ ] = 0, it is clear 

that (Hi - E^)^ = as required. The SUSY charge 
operators are also vectors with components defined by: 

Qi = di + W i: Q\ = -d t + W t . (II.7) 

We can now define the sector two Hamiltonian such 
that, above the ground-state (Eq), it is isospectral with 
Hi. We do this as follows, for the first excited state in 
sector one we can write 

QlQ^ = (E[ 1) -E { 1) )4 1) . (11.8) 

We then form the tensor product by operating on the left 
with Qi so that 

(QQt) . Q^W = (£« _ E^)Q X ^\ (II.9) 

That is to say using Einstein notation, 

(QiQ^Qj^ = (e[ 1] - 4 1} )Qi4 1] - ( n -l0) 

It then follows that Qipi^ is an eigenstate of the tensor 
Hamiltonian W2 = (QQ^) with energy E^ = e[^ — 
Eq . Since we are free to set the energy origin, taking 
E { 1} = gives E { 2) = E^\ It is also clear that Qtp^ 
cannot generate a lower energy eigenstate of ^2 since 



3 



Qib^ = 0, so that Qipi is indeed proportional to the 

ground state of n 2- The precise relation between the 
two sets of states is given by 
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(i) 

n+l 



(i) 



III. SUSY-QM FOR THE 
THREE-DIMENSIONAL HYDROGEN ATOM 

We now consider the hydrogen atom. We begin by 
noting that the ground state is exactly given by 



^1,0,0 



(III.l) 



where we have set the Bohr radius equal to 1 . The Hamil- 
tonian (in atomic units) is simply 



H 



1. 



Then, the vector superpotential is given by 
W = -V In -01,0,0 = f, 



(III.2) 



(III.3) 



where f is a unit vector in the direction of f. This is 
an extremely interesting result. First, we see that the 
superpotential for the Coulomb interaction is, itself, non- 
singular. Second, in the standard approach, because the 
angular degrees of freedom have already been separated 
out, the superpotential is a scalar and it depends on the 
angular momentum squared (i.e. on 1(1+1)). The precise 
form for the ground state (Z = 0) is 



and 



v 2 mi 



(11.11) yields 



V 2 0i, o ,o 



2E Q 



"01,0,0- 



(III.9) 



(III. 10) 



Since the ground state energy of hydrogen in atomic units 
is -1/2, we find that Equations (III. 7 1 and (III. 10 1 are 



consistent and we have obtained the correct vector su- 
perpotential. Of great interest is the wave equation for 
the sector two problem. This Hamiltonian is given by 



«-> 1 

n 2 = -^vv 



WW + VW 



(III. 11) 



In the case of the hydrogen atom, because we have exact 
analytical expressions for the excited states of Hi, it is 
a simple matter to generate analytical expressions for all 
the states of the sector two Hamiltonian. It is convenient 
to label the sector two states with an index indicating 
the n th energy state (i.e., we use the principle quantum 
number n = 1,2,...) along with the quantum numbers of 
the sector one excited state from which they are obtained. 
Thus, the four degenerate ground states of will be 



denoted by Vi.ap*. V4,iv Wx,L- We choose here 

to use the real states rather than those labeled by mi — 
±1 and mi = values. We find that these solutions are 
given by 



^ (2) 



N 



-r/2 



-r/2 



(III. 12) 



W ra dial — 1 

In three dimensions, we have 

W = el- + €:,- +e. 



(III.4) 



(III.5) 



The magnitude of W is equal to the radial superpoten- 
tial, as one expects, but the individual components are 
radically different. Note that these components can also 
be written solely in terms of angular functions (the di- 
rection cosines of f). To obtain the atomic potential for 
hydrogen, we form 



w -W -V -W = 1 



V 



-2E Q 



2 



(III.6) 



(III.7) 



$ A* = N 



N 



k e -r/2 + Z J_ e - 



r/2 



4% = -^2 e ~ r/2 - 



(III. 13) 



(III. 14) 



(III. 15) 



These equations can be verified by simply applying Q 
to the first excited state wave functions of sector one. 
It is also easily verified that acting on these states 
regenerates the ip^ and ip^J states. Furthermore, in 
Figures 1 and 2 we provide plots of the i^ls an( ^ Yi\p x ' 

7*(2l ~^(2) 

It is straight forward to see that 1^ and ip{ 2 Pz are 
both similar to 0i,2 Px - 



Now we recall that 



H4> mi = E n ip mi 



(III.8) 
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FIG. 1. The three components of the wave function for ip\ s- Here, blue corresponds to positive values and red to negative. 



I Pi Component. I] I Pi Component. 2] [Pi Component, 3] 




-*(2) 

FIG. 2. The three components of the wave function for ip\ ' Here, blue corresponds to positive values and red to negative. 



5 



IV. AN APPROXIMATE SUPERPOTENTIAL 
FOR THE HELIUM ATOM 



V. AUFBAU APPROACH FOR EXCITED 
STATES 



It is of interest to begin exploring how our approach 
to multidimensional SUSY-QM would deal with a two 
electron atom. It is clear that the usual radial (one- 
dimensional) hydrogen atom SUSY-QM treatment is not 
readily generalizable to deal with helium. We have car- 
ried out a Quantum Monte Carlo study of the sector one 
ground state of helium using the Pade Jastrow trial wave 
function: 



$1 = e-^e-^t 



2(1 + ° 



(IV.l) 



with the optimum a given by a = 0.353. This yields 
an energy of e[ 1] w 2.878, which is in error by about 
1%. This error is reasonable for a simple treatment ne- 
glecting relativistic interactions. The approximate W is 
generated from 



-2ri - 2r 2 



ri2 



2(1 + ar 12 ) 



(IV.2) 



(IV.3) 



Here 
V = 



f lx 



d 



„ d „ d d _ <9 „ d 

Ely Tj r-eis TJ \-&2x Tj \~^2y Tj l~ e 2z 7i > 

dt/i dzi cte 2 cty 2 Cz 2 

(IV-4) 

where the {ey} are orthonormal vectors. The resulting 
vector superpotential for the Pade- Jastrow trial function 
is readily found to be 



W He{PJ ) = 2fi + 2f 2 - fia 



1 - 



(1 + ar i2 ) 



(IV.5) 



Thus, the structure of Wh&{p,j) is analogous to Wh in 
that Coulomb interactions generate vector superpoten- 
tials that involve unit vectors anti-parallel to the direc- 
tion of the forces. This is true in general for Coulom- 
bic interactions. This emphasizes the important dis- 
tinction between our three-dimensional SUSY-QM treat- 
ment of an atom and the standard hydrogen atom one- 
dimensional radial SUSY-QM. We are currently carrying 
out calculations of the sector two ground state for he- 
lium using the above Wue(pj)- There are several possi- 
ble approaches to be explored. One is to follow our ear- 
lier two-dimensional study and employ a Rayleigh-Ritz 
variational method[l . This is immediately applicable to 
the helium problem and should work without difficulty. 
However, it will involve much more computational effort 
since helium is a six-dimensional system as opposed to 
the two-dimensional systems studied earlier. The sec- 
ond approach we intend to explore is the Dirac-Frankel- 
McLachlan time-dependent variational method|22]. The 
results of these studies will be reported later. 



For multielectron atoms, it becomes necessary to con- 
sider how the aufbau principle acts in the second sector to 
permit efficient calculations of sector one excited states. 
This is because we can use this principle to design rea- 
sonable trial wave functions for a variational approach to 
the sector two ground state. In this section, we consider 
a simple aufbau description of the sector one helium ex- 
cited states in order to design an approximate sector two 
ground state of helium. We assume that in the first ex- 
cited state of sector one, we have one electron in the 15 
orbital, given by wave function a, and one electron in the 
2S orbital, given by wave function /3 where 



-2r 



a(r) 



and 



4V2tt 



(1-r). 



(V.l) 



(V.2) 



Then, it is of interest to take the product of these states 
such that we have a(ri)/3(r 2 ), to which we can apply our 
Q to find 

Q(a(n)/3(r 2 )) = - e - 2r ^ [2h (1 -r 2 ) + f 2 ] = <^ 2) . 

(V.3) 

It is clear that to find the second sector state associated 
with a(r 2 )/3(ri), we simply need to interchange labels 1 
and 2 in Equation (V.3 1 to get 



Q(a(r 2 )/3(n)) 



e -2r 2 - ri [2f 2 (l-n)+f 2 ] = P 12 <^\ 

(V.l) 

where P 12 exchanges the electron labels. Then, we can 
use this as a "building block" to construct our ground 
state in the second sector in the following way. Suppose 
we want to create the the spin triplet state, we find that 
we need to subtract the first building block from the sec- 
ond. This gives us a second sector result of 



-1r\ —T2 



[2fi (1 - r 2 ) 



T2\ 



-2r 2 — n 



pfaCl-rO+fi] 



(V.5) 



And similarly, we can find the second sector electronic 
structure singlet state by simply adding the two building 
blocks as given below. 



$?L ? <et=-e- 2n - r2 [2fi(l-r a )+f 2 ] 

£ -2r 2 - ri [ 2f2(1 _ ri)+fl ]_ 



(V.6) 



Indeed, by taking the scalar product with Q\ we can 
verify that ipfj ripl?t and ^ } singlet give the appropriate 
spatial wave functions. This is to say that, to within a 
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multiplicative constant, we get that 

Q f ' triplet = # = <n)p{r 2 ) - a{r 2 )p{r x ) (V.7) 
and 

Q f ■ i?L g iet = 4 1} = «(ri)/?(r 2 ) + a(r 2 )/3(n). (V.8) 

From this, we observe that the aufbau principle in the 
second sector is remarkably simple. We merely need to 
take the building block <p\ and antisymmetrize or sym- 
metrize appropriately. One interesting thing to note is 

(2) 

that our "building block", cp\ , is neither symmetric nor 
antisymmetric under particle exchange. 

However, this basis doesn't include the correlation. To 
do this, we can multiply our antisymmetrized second sec- 
tor wave function by a correlation function, given by the 
Pade-Jastrow function which only depends on ryi- It is 
clear, then, that because our correlation function is only 
a function of ri2, its symmetry will not be affected by 
the application of Q and, thus, we can simply multiply it 
by our second sector state of interest (where the minus 
corresponds to the triplet and the plus to the singlet): 

triplet =e^k*r (-e- 2r ^ [2r x (1 - r 2 ) + r 2 ] T 
e -2r 2 - ri [2f 2 (l_n)+ri]). 

(V.9) 

These equations appear to be qualitatively correct but 
one will insert variational parameters (e.g., effective 
charges, etc.) when doing computations. The next step 

is to perform numerical calculations using ipl l rip i et in the 
second sector to find energies generated only from the as- 
sumptions above and compare the results with standard 
approaches. We are currently exploring more complex 
atomic structures to generalize this further. 



VI. CONCLUSIONS 

In this paper we have shown how our multi-dimensional 
generalization of SUSY-QM can be applied to the hydro- 
gen atom. Previously, most detailed attempts to treat 
the hydrogen atom first separated the angular degrees of 
freedom, leaving a one-dimensional radial wave equation. 
It was then possible to obtain the SUSY-QM factoriza- 
tion, yielding a scalar superpotential that, for the / = 
states, is simply W — 1. While these results are inter- 
esting, the one-dimensional radial SUSY-QM approach is 
not readily generalizable to treat even the helium atom. 

In our approach, the full three-dimensional character 
of the hydrogen atom is considered, with the result being 
a vector-valued superpotential, W , which for the hydro- 
gen atom, is W = r. That is, the vector superpotential 
points in the opposite direction of the attractive Coulomb 
force between the electron and the nucleus. This is inter- 



esting also because, although the Coulomb potential is 
singular, its vector superpotential is not. It is important 
to note that such a superpotential was also obtained ear- 
lier by Stedman [2j. However, his sector two Hamiltonian 
differs from ours and produces "extra" states that are not 
degenerate with sector one. 

The fact that W for the three-dimensional hydrogen 
atom is a vector does not, in any way, modify the sec- 
tor one dynamical equation. However, the sector two 
situation is radically affected! In the one-dimensional 
SUSY-QM case, there is no significant change in the 
basic mathematical structure of the sector two partner 
Hamiltonian. In the multi-dimensional case, the sector 
two Hamiltonian is a tensor. However, we have shown 
in previous studies, that many of the standard compu- 
tational techniques remain valid. Of particular inter- 
est is the Dirac-Frankel-McLachlan Variational Method, 
since this is known to deal better with higher-dimensional 
systems |22|. 

In the case of the hydrogen atom, it is straight forward 
to generate all the sector two eigenstates. This is a con- 
sequence of the fact that exact analytical eigenstates of 
the three-dimensional hydrogen atom are known. It is 
then easy to apply the charge operator, Q, to the excited 
hydrogen atom states and obtain sector two eigenstates. 
(We note that because of the four fold degeneracy for 
the sector two ground state, the resulting eigenstates can 
be super-posed in any manner convenient for the study 
at hand). It is of considerable interest to begin explor- 
ing how our multi-dimensional SUSY-QM treatment can 
be applied to the helium atom. In this case, the ex- 
act sector one ground state is, of course, unavailable. In 
our previous one and two-dimensional studies we have 
considered other systems for which an exact W was not 
possible. In the case of helium, we chose to examine an 
accurate Pade-Jastrow approximation to the sector one 
ground state. In this case, it is easy to obtain an analyti- 
cal (albeit approximate) W that displays very reasonable 
intuitive character. In direct analogy with the exact hy- 
drogen atom W, we find that the Wne(pj) vector super- 
potential consists of a combination of unit vectors that 
again, are anti-parallel to the Coulomb forces associated 
with the helium atom potential energy. The next step 
in our study will consist of computations of a sector two 
ground state, which will allow us to obtain an approx- 
imate helium atom sector one first excited state energy 
and wave function. 

Future studies will explore extending the approach to 
more than two electron atoms. There, the issue will be 
taking account of the electrons' spin degrees of freedom. 
Our current plan is to employ the "spin-free" techniques 
of Matsen p3H55| and others |51). 

We have also generalized the aufbau principle to work 
in the second sector Hamiltonian, demonstrating that we 
are able to produce reasonable forms of excited states 
by simply using hydrogenic orbit als. The equations have 
a reasonable structure but variational computations are 
necessary. We shall report these results later. 
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